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ABSTRACT: General properties of new models of the electronic Fock spaces 
based on the notion of {p, g)-sheaves are studied. Interrelation between simple 
sheaves and density operators is established. Explicit expressions for the trans- 
formed reduced Hamiltonians in terms of the standard creation-annihilation op- 
erators are presented. General scheme of parametrization of p-electron states by 
xr-electron means (>f = 2,3,...) is described and studied in detail for the case 
of sheaves induced by >f-electron wavefunctions. It is demonstrated that under 
certain conditions p-electron problem may be reformulated as the eigenvalue prob- 
lem in >f-electron space equipped with certain p-electron metric. Simple numerical 
examples are given to illustrate our approach. 
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Introduction 

In present work we continue study of properties of new models of p- 
electron sections of the Fock spaces introduced in our previous work 
These models are based on the notion of {p, g)-sheaves that can be considered 
as g-electron representation of p-electron states [q < p). Figuratively speak- 
ing, the space of all (p, g)-sheaves is a model of p-electron space 'ready for 
g-electron interactions'. This means that no contraction is required for calcu- 
lation of matrix elements of arbitrary g-electron operators. The contraction 
is replaced by the summation over g-electron functions (germs) constituting 
the sheaf. 

In the second section necessary basic definitions are given. 

In the third section general properties of simple {p, g)-sheaves are studied. 
Important notion of vector Z-ceW is introduced. Explicit characterization 
of imagies of simple sheaves with respect to different compositions of the 
assembling-disassembling mappings is given. 

In the fourth section it is shown that there exists a one-to-one correspon- 
dence (up to arbitrary overall phase prefactor) between simple sheaves and 
a certain class of density operators. 
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In the fifth section the explicit expression for the transformed reduced 
Hamiltonian in terms of the creation-annihilation operators is obtained. 

In the sixth section the notion of p-electron metric in >i:-electron space 
(xr = 1,2,3,.. .) is introduced and it is demonstrated that under certain 
conditions p-electron problem reduces to the eigenvalue problem in >ir-electron 
space equipped with p-electron metric. 

In the seventh section orbital representation of (p, g)-sheaves is discussed 
and relevant formulas for metrics and Hamiltonians are derived in terms of 
molecular orbitals. 

The eighth section is dedicated to numerical testings. 



Basic Definitions. 

For fixed basis set of n orthonormal spin orbitals the corresponding finite- 
dimensional Fock space Tn is spanned by determinants \R) where R runs over 
all subsets of the spin-orbital index set. Basis determinants will be labelled 
by subsets and all sign conventions connected with their representation as the 
Grassman product of ordered spin-orbitals will be included in the definition 
of the creation-annihilation operators and in the definiton of specific set- 
theoretical operation Ak, K G N that was introduced in and studied in 
detail in 0.Here is the spin-orbital index set. 

For any {p + g)-element subset Z of the index set let us denote by the 
symbol the following subspace of the g-electron sector of the Fock 

space: 

(9) 

^n,,(^) = 0C|S) (1) 

scz 

where C is the field of complex numbers. Let us define the set 

Bn,p,g = {{Z, S) C N X N : \Z\ = p + qk\S\ = qLS C Z}, (2) 

and the equivalence relation on this set 

{Z, S) ~ {Z', S') ^ Z\S = Z'\S'. (3) 

The set of the corresponding equivalence classes contains elements and 
in each equivalence class there are elements. 
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Definition. Family {ipz}zcN of g-electron functions, ipz ^ ^n,q{Z), is 
called a sheaf of q-electron germs of p-electron wavefunction, or just a (p, q)- 
sheaf, if the mapping 

{Z, S) ^ (-l)l(A5)nA.| ^^1^^^^ \Z\=p + q, \S\ = q,ScZ (4) 

is constant on the equivalence classes of the set Bn,p,q modulo the equivalence 
relation (3). Here Az is the mentioned above set-theoretical operation that 
is used to reduce manipulations with sign prefactors in Grassman algebra to 
pure set-theoretical ones. 

The set of all (p, g)-sheaves is denoted as Sn,p,q- 

In paper it was demonstrated that there exists a one-to-one correspon- 
dence between the set of all (p, g)-sheaves and vectors of the p-electron sector 
of the Fock space. The mapping 

(p) 

^n,p,q ■■ Cr\R) {^z}zcN (5) 

RCN 

where 

(g) 

\^z) = 5^(-l)l(^\^)"^-IC7^\5|5). (6) 

scz 

is called the disassembling mapping and is used to transfer all relevant struc- 
tures from the p-electron sector of the Fock space on the set of all {p, q)- 
sheaves. The linear structure is transfered as 

5„,p,,($^A.vI;«) = =J2^^i'^z}zcN. (7) 

i Li J ZcN i 

This equality just expresses the fact the due to the linear character of the 
gluing conditions (4) arbitrary linear combination of {p, g)-sheaves is a {p, q)- 
sheaf. 

The inner product (Euclidean structure) is defined as 

{p+g) 

{{4^z}zcn\{<1^z}zcn) = Yl ^^^l*^^) 

\ q ) ZcN 

and it is consistent with the inner product of the corresponding p-electron 
wavef unctions. 
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The inverse to Sn,p,q is called the assembling mapping and is used to restore 
|}-electron vectors from (p, g)-sheafs with the aid of the gluing conditions (4). 

In paper it was mentioned that from the very beginning it seems rea- 
sonable to associate with each p-electron wavefunction \l/ its sheaf Sn.p.q{^) 
and co-sheaf Sn,n-p,q{p^) where I : \R) — > |A^\-R) is the particle-hole in- 
volution. However, since general properties of sheaves and co-sheaves are 
identical, and actual difference appears only on the level of evaluation of re- 
duced Hamiltonian matrix elements, we study properties of sheaves, keeping 
in mind that the analogous properties are valid also for co-sheaves. 

The family 

{7r„,p,g(Z) : {iJz'}z'cN i'z}zcN (9) 

involves linear mappings that perform projections of the vector space Sn^p^q 
on the vector spaces Tn,q{Z\ 

Let us suppose that some g-electron wavefunction '^z ^ ^n,q{Z) is chosen: 

(g) (p) 

i;z = J2^s\S) = (10) 

scz RcZ 

where R = Z\S and 

Cn = (-l)l^^^-lcs (11) 

Since J-'n,q{Z) is isomorphic to Tn,p{Z\ we can easily lift g-electron wave- 
function to the level p using, for example, the linear mapping: 

(p) _ (P) 

u\{Z) : 5^(-l)l«^^-IC«|Z\i?) ^Y.^r\R) (12) 

RcZ RcZ 

Subsequent disassembling of p-electron wavefunction vP^{Z^'^z gives us 
(p, g)-sheaf that may be considered as generated by g-electron wavefunction 
'^z- Compositions Sn^p^q o u^{Z) = in,p,q{Z) constitute the family of linear 
injective mappings 

in,p,q{Z) : T n,q{Z) Sn,p,q (13) 

that are right inverses (sections) of the projections (9). 

Definition, (p, g)-sheaf generated by g-electron wavefunction il)z is called 
simple and is denoted as {iIjzz'}z'cn- 

It is to be noted that the same simple (p, g)-sheaf may be generated by 
different g-electron wavefunctions. 
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{p, g)-sheaves corresponding to single determinant p-electron wavefunc- 
tions are called determinant sheaves. The characteristic property of deter- 
minant sheaves is: simple {p, g)-sheaf is the determinant one if and only if it 
is generated by any of its nonzero germ (see|2]). 

It is clear that the assembling and disassembling mappings may be used 
to establish isomorphism between operator spaces over J~n^p and over S^^p^q. 
Indeed, if /i G Endc{Tn,p) then the mapping 

h Sn,p,q oho s'l^g (14) 

is the required isomorphism. For physically relevant operators it is usually 
possible to express their matrix elements between two {p, g)-sheaves as a sum 
of matrix elements of certain g-electron operators between germs of these 
sheaves. In our previous paper the electronic Hamiltonian H was handled 
in such a way. Namely, due to the specific form of H, its matrix element 
between two {p, g)-sheaves may be written as 

. ip+g) 

{{4>z}zcn\H\{iPz}zcn) = rn-v\ 

\ q ) ZCN 

where 



p _ 1 2 ri2 



(16) 



is the so-called reduced Hamiltonian and A{n,p, q) is introduced in oper- 
ator acting on the operator space JF„ ^ (g) JF* ^. 



Properties of Simple Sheaves. 

With each function ipz £ ^n,q{.Z) written in the form of Eq.(lO) let us 
associate two sets: 

/p(V>z) = {R^Z:\R\=pLCr^^] (17) 
L{i^z)= U R (18) 

Let {'ipzz'}z'(zN be some simple (p, g)-sheaf generated by g-electron wave- 
function ipz ■ Then any germ of this sheaf may be written as 
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te'= E (19) 

R'<zz'r\z 

where the coefficients Cri are given by Eq.(ll). 

Proposition 1. ^-electron wavefunction generating simple (p, gf)-sheaf is 
necessarily a germ of this sheaf. 

Proof. Follows directly from the relation 

where idjr^^^(^z) is the identity operator over Tn,q{Z)M 

Proposition 2. g-electron functions '^l)z-^ and -0^2 generate the same 
simple (p, g) -sheaf if and only if 

(i) Wz^) = Wz^); 

(ii) For any it! e Ip{ipZi) functions ipzi and ipz^ are subject to the following 
"gluing" conditions 

Proof. Let us write down g-electron functions ipzi and ipz2 in the form 
of Eq.(lO) 

(p) 

V'z. = Yl (-1)'^"^^^I4^^IM^) = 1,2) 

RcZi 

From the equality {V^ZiZ'jz'civ = {ipZ2Z'}z'cN it follows that ipzi = ipz^Zi 
and ■(/'^j = ipZiZ2- These equalities together with Eq.(19) imply that C^jp — 
and, consequently, condition (i) and (ii) are fulfilled. 
Now let us suppose that for given functions ipzi and ■0^2 the conditions 
(i) and (ii) hold true. Then from (i) it follows that for arbitrary Z' C N 

(p) 

^z,z'^ E (-i)'^"^-''4^m^), 

RcZ'nZinZ2 

where i=l,2. Use of conditions (ii) immediately leads to the equality 

{'ipZiZ'}z'cN = {4'Z2Z'}z'cN ■ 

Proposition 3. The number of g-germs generating simple (p, q)- sheaf 
{i/^zz'jz'cN is equal to 

/ n-\L{^z)\ \ 
\p + q-\L{tl;z)\J 
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Proof. By definition, L{%l)z) C Z . For any Zi D L^ipz) the germs ijjz 
and ipzzi obviously satisfy the conditions of proposition 2 ■ 

Corollary 1. For any determinant sheaf |L('j/^2)| = p and it is generated 
by any of its {"'~^) germs. 

CorollEiry 2. Non-determinant (p, l)-sheaf has only one generator. 

Proof. For non-determinant sheaves \L{il)z)\ should necessarily be equal 
to p + 1 ■ 

Definition, g-electron functions ibzi and ibz2 are called p-equivalent 
{i^Zi ~ "^22) if thsy generate the same simple (p, g)-sheaf. 

Proposition 4. Let {'^zz>}z'<:in be simple (p, g')-sheaf generated by q- 
electron wavefunction tpz- Then any g-electron function from Tn,q{Z) is 
orthogonal to all ipzz' with Z' 7^ Z. 

Proof. Expansion of arbitrary g-electron function from Tn.q{Z) (see 
Eq.(lO)) involves g-electron determinants \Z\R) where R <Z Z. From Eq.(19) 
it follows that each g-germ i^izz' of the sheaf under consideration is expanded 
via g-clcctron determinants \Z'\R'). It is easy to see that Z' ^ Z and 
R' d Z'r\Z implies Z'\R' gL Z m 

Corollary, g-germs generating simple (p, g)-sheaf are mutually orthogo- 
nal. 

Let us consider a family of k simple (p, g)-sheaves {ipz[z'}z'(iN-i i — 
1,2, ... ,k and sort out this family in the following manner: 

Select some subset Zi from the family Z[, Z'2, . . . ,Z'^ and collect all (p, q)- 
sheaves such that "ipz^Zi ~ V'zn The resulting subfamily will be denoted as 

{'^zlz'}z'cN,h = 1,2, 

From the rest of the initial family select subset Z2 and collect all (p, q)- 
sheaves such that ipz'.Z2 ~ The resulting subfamily will be denoted as 

{i^zlz'}z'cN,i2 = l,2,...,k2. 

Repeating the last step till the initial family is exausted, we arrive finally 
at the family {ip^ZiZ'^z'cNiii = 1,2, . . . , ki,l = 1,2,...,. By construction, 
each subfamily of this family corresponding to some fixed index 1 involves 
simple (p, g)-sheaves generated by g-germs from the vector space Tn,q{Zi). 
As a result, the problem of orthogonalization of simple sheaves within this 
subfamily seems to be not very complicated, especially if one takes into ac- 
count the following assertion. 

Proposition 5. Let {tpzz'}z'cN and {'^^I'I-Z'cat be simple (p, g)-sheaves 
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generated by ^-electron functions ijjz\'^^z^ £ ^n,q(Z). Then 



(20) 



Proof. Direct calculation based on the definition (8) of the inner product 
leads readily to the equality required ■ 

The question when linear combination of simple sheaves is simple is of 
a certain interest in connection with the problem of selection of linearly 
independent simple sheaves. The answer to this question is given by the 
following assertion. 

Proposition 6. Linear combination 



i=l 



of simple (p, g)-sheaves is simple if and only if there exists (p + g)-element 
subset Z C N such that 



where 



TZl...Zk 



{l<i<k: ZiD R}. 



Proof. Arbitrary linear combination of simple sheaves may be written 
in the form 



(p) 



/ 



1=1 



R'dZ' 



Z'CN 



If we suppose that this sheaf is simple then there should exist index set 
Z G N such that the g-germ 



RCZ \^e4-'>' J 



generates this sheaf. It can, however, be true if and only if the coefficients 
by the basis determinants \Z'\R') in expansion of gr-germs are equal to zero 
for all R' <;tZU 
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Corollary. Arbitrary linear combination of determinant {p, g)-slieaves 
generated by g-germs 

k 

is simple if and only if there exists index set Z C iV such that ZD IJ -^j) 

k 

or, in other words, | IJ -Rj] <p + q- 
1=1 

Proof. Without loss of generality we may assume that the coefficients Aj 
of linear combination of determinant sheaves are non-zero. Let us suppose 
that for any subset Z d N there exists Rj (f_ Z. In this case 



Zl...Zj. 



and, consequently, the linear combination under consideration can not be 
simple ■ 

Definition. The image of JF„ ^(Z) with respect to the monomorphism 
jn,p,q{Z) is called the vector Z-cell of the vector space of all {p, g)-sheaves 
and is denoted as Jn,p,q{Z). 

It is pertinent to note that Z-cells depend on the MSO basis chosen. 

Each vector Z-cell is a (^^"^j -dimensional subspace of Sn,p,q- The set of 
all simple (p, g)-sheaves in Sn,p,q is a union of its vector Z-cells: 

(p+<?) 

Jn,p,q ~ \^ •Jn,p,qiZ) . (21) 
ZCN 

It is easy to show that the set Jn,p,q generates the vector space of all (p, q)- 
sheaves (see |]T|). 

The set of simple (p, g)-sheaves {4'zz'}z'cn ^ Jn,p,q{Z) such that 

\L{^z)\<k 

where k = p + q,p + q — 1, . . . ,p, will be denoted as Jn},^q{Z) . It is clear that 
each Z-cell admits the following filtration 

JnUz) = Jiti^z) ^ J^i^r'^z) ^ • • • ^ J^nU^) (22) 
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Topological dimension of Jn%q(Z) is equal to (^). Thus, each Z-cell possesses 
fc-faces (borders) where simple (p, g)-sheaves from different cells are situated. 
It is likely that filtration (21) may serve as a starting point for further more 
advanced analysis of geometry of the set Jn,p,q- 

Now let us consider the connection between spaces Sn^p^q with different q. 
There exist obvious isomorphisms 

^n,p,q ^ n,p ^ ^n,p,g±l- v / 

Proposition 7. The image of any simple (p, g)-sheaf from Jn,p,q with 
respect to the isomorphism Sn,p,q+i o s^j, ^ is simple (p, q + l)-sheaf, and for 
k ^p + q,p + q-l,...,p 

-'^n,p,q+l ° ^n,p,q {'^n,p,q) ^ '^n,p,q+li i'^^) 

that is the image of Jn,p,q belongs to the border of Jn.p.q+i- 

Proof. Let ipz ^ ^n,q{Z) generate simple {p,q)- sheaf {ipzz'jz'cN and 

(p) 

u^{Z)il)z — Cr\R) — "ifz- For any index i e N\Z the p-electron wave- 

RCZ 

function 

(p) 

*zu{i} = CR\R) = 'ifz 

RCZUi 
(i?3i) 

can be disassembled to give simple (p, q + l)-sheaf with germs 

(p) 

i^zu{i}z'= Yl (-i)i''"^^''c'i^|z'\i?) 

Rcz'nz 

where \Z'\ = p + q + 1. It is clear as well that the choice of concrete index 
i is irrelevant. Indeed, for any index j e N\Z, {q + l)-germ ipzuiijzuij} is 
p-equivalent to ipzuii} ■ 

Lowering of simple (p, 5)-sheaf to the level q — 1 leads, in general, to linear 
combinations of simple (p, q — l)-sheaves. 

Proposition 8. Let {ipzz'}z'cN be simple (p, g)-sheaf generated by q- 
germ tjjz- Then 

Sn,p,q-1 O Sn,p,5 ({fc'}z'Civ) = " ^{V'Z\{fc}Z"}z"CiV- (25) 
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Proof. Disassembling of p-electron function — u^iZ)ijjz, \Z\ — p + q, 
gives (p, q — l)-sheaf with germs (see Eq.(6)) 

(p) 

i^z"= Yl {-iy'^'^'"^CR\Z"\R), (26) 
Rcz"nz 

where \Z"\ = p + q — 1. 

On the other hand, let us consider the sum 

Y{'4^z\{k}z"}z"cN = I ^ipzMkjz" > , (27) 
fcez Kkez ) z"cn 

where simple (p, q — l)-sheaf {'4'z\{k}z"}z"cN involves {q — l)-germs 

(p) 

MWZ"= E {-iy^^^^"^Cn\Z"\R). 

Rc[Z"niz\{k})] 

Thus, (p, q — l)-sheaf on the right-hand side of Eq.(27) involves {q — 1)- 
germs 

(p) 

J2 Yl (-i)i''"^^"ic'fl|z"\i?). 

kez Rc[Z"n{z\{k})] 

The contribution to this sum from k ^ Z"r\Z is equal to {p + q— \Z" r\Z\)ipz" 
where ipz" is given by Eq.(26). On the other hand, for each fixed R C 
Z" n {Z\{k}) there exist \Z" D Z\ - p different k E Z" f] Z. As a result, 

^{V'z\{fe}Z"}z"cAr = q {tpz"}z"cN ■ 

fcez 

It is to be noted that if L{iIjz) is a proper subset of Z then for any 
k e Z\L{tpz) (the choice of k is irrelevant), Sn,p,q-i z) = Sn,p,q-i z\{k]) 
is simple ip,q — 1) -sheaf. 

Corolleiry. For any admissible q and any x < q 

«n,p,^ O Sn!p,q {{tpZZ'} Z'cn) = 3 {'^Z\KZ"} Z"cN ■ (28) 

^' KCZ 

Very important is the following statement that is contrary in a certain sense 
to the Proposition 8. 
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Proposition 9. If Z C N, \Z\ — p + q, and for each {q — x)- element 
subset K <ZZ 

(p) 

i>z\K= Yl {-iy'^'^''^''^C^R^\Z\K\R), (29) 

RCZ\K 

then the image of the sum 

{'^Z\KZ'}z'<zN 

K<ZZ 

with respect to the assembhng mapping s^j, is p-electron wavefunction 

(P) / (9--) \ 

= E E (30) 

R<zz \kcz\r J 

and its disassembling with the aid of Sn,p,q gives simple (p, g)-sheaf generated 
by the ^-germ 

(p) / \ 

^(_l)|HnA.| ^ (31) 

Proof. By direct calculation. 

This Proposition can be used to parametrize subsets of the sets Jn,p,q{Z) 
using x-electron means (x < q). 



Simple Sheaves and Density Operators 



Definition. The quadratic mapping dq : Sn,p,q — > J^n,q ® ^n,q defined by 

. (p+q) 

dq {{i^z}zcN) = 7^ E (32) 
\ g ) ZcN 

is called the density mapping of order q. The value of this mapping at some 
(p, g)-sheaf is called the g'-density operator associated with this sheaf. 

The problem of description of the set dq {Sn,p,q) is called the pure repre- 
sentability problem and in slightly different form was formulated in the very 
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first papers by Coleman ^ ^, |^. By an abuse of notation tlie same symbol 
dq will be used for the density mapping defined by Eq.(32) and for elements 

of the set dq {Sn,p,q) 

Let us consider density operators corresponding to simple sheaves. From 
Eq.(19) it follows that 



Z'CN) 



\i^z){i^z\ + 



ip+q) 

E 



\2p{q) 
I ^N\Z 



Z'cN 

(p<\z'nz\<p+q) 



(33) 



where 



P 



(<?) 



X 



(34) 



(g) 
Sex 

is the projector associated with the index set X G N. Note that the rep- 
resentation (33) of the density operator associated with some given simple 
sheaf is not unique (see Proposition 3). 

Complete characterization of g-density operators corresponding to simple 
{p, g)-sheaves is given by the assertion that readily follows from Eq.(33) and 
Proposition 4. 

Proposition 10. g-electron operator dq is associated with simple {p, q)- 
sheaf if and only if there exists {p + g)-element subset Z C N such that 

(i) the operator P^^dqP^^ corresponds to some (not necessarily normal- 
ized) pure g-electron state, that is P^^dqP^^ = \ipz){ipz\', 

(ii) (/ - ) dq (l - P^^^) = 

V /V / Z'CN 

(p<\z'nz\<p+q) 

where wavefunctions ipzz' are uniquely determined by the eigenfunction 
ipz of the operator P^^dqP^^ corresponding to its non-zero eigenvalue. 

Thus, for any fixed {p + g)-element subset Z <Z N there exists one-to-one 
correspondence (up to arbitrary overall phase prefactor) between elements 
from the vector Z-cell Jn,p,q{Z) and representable density operators of the 
form of Eq.(33). Reservation 'up to arbitrary overall phase prefactor' may 
be removed by turning to the quotient of Jn,p,q{Z) obtained by identification 
of sheaves that differ by overall phase prefactor, or to the projective Z-cells 
if only normalized sheaves are of interest. We, however, prefer to use vector 
Z-cells saying that simple sheaves parametrize density operators of the form 
of Eq.(33) even if this statement is not quite correct. 
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The simplest density mapping corresponds to the case q — 1 and its value 
at some simple {p, l)-sheaf is equal to 



n — p 



n — n ii \^ 



(35) 



where \Z\ = p + 1 . The kernal of the operator on the right-hand side of this 
equation is p-dimensional subspace of occupied one-electron particle states 
(occupied molecular spin orbitals). 

Using proposition 7 and Eq.(33), it is possible to lift one-density operators 
from di[Jn,p,i) to the level q = 2: 



d2{{'ipZlJ{i}Z'}z'cN) 



C7) 



|V'ZU{j})('0ZU{i}l 



(P+2) 



+ 



J2 I ^ZU{i}Z' ) {ipzu{i}Z' I + Pn\^zu{ 



i}) 



Z'CN 

(|z'nz|=p+i) 



(36) 



Non-zero contributions to the sum on the right-hand side of Eq.(36) may 
appear in two cases: Z' = ZU {j},j G N\Z,j ^ i and Z' = Z\kU {i,j},k G 
Z, i,j G N\Z. Analysis of these cases leads to the following final expression: 



d2{{lpZlJ{i}Z'}z'cN) 



where 



1 



C7) 



k,k'ez jeN\z 



(2) 
N\Z 



(37) 



0, if A; = k' or k, k' > j; 
e = I [{Z\{k})A{Z\{k'})]nA^,}\ = {1, iik<jkk'>jork'<jkk> j; 

2, if k, k' < j. 

By the symbol {j,k} on the right-hand side of Eq.(37) non-ordered pairs 
of spin-orbital indices (two-element subsets) are denoted. It is pertinent to 
emphasize once again that the right-hand side of Eq.(37) does not depend 
on the choice of index i G N\Z that appears explicitely on the left-hand side 
of this equality (see Proposition 7) . 
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It is possible, following Coleman, to replace the requirement of pure repre- 
sentability by more simple requirement of representability by an ensemble of 
pure states. Instead of the quadratic g-density mapping (32) let us introduce 
the mapping 

that is defined on generators \{4'z}zcn){{<Pz}zcn\ as 

\^z){<Pz\ (38) 

\ q ) ZCN 

and is continued to Sn,p,q ® ^n,p,q by linearity. 

The set of all finite convex combinations of the type 

ZCN\ 

i 

where Aj > 0, ^ Aj = 1, and | KV'^^jzcAfl I = 1 will be denoted as Ens{Sn,p,q)- 

i 

Its image with respect to dq is the set of all ensemble representable g-density 
operators and the problem of its analytic description is known as the ensemble 
representability problem. For g = 1 its constructive solution was found by 
Coleman [|[. 



Explicit Expressions for the Transformed Reduced Hamiltonians 



From Eq.(15) it follows that in the vector space of all (p, g)-sheaves the 
usual p-electron Hamiltonian should be replaced by the reduced Hamiltonian 
transformed with the aid of the automorphism A~^{n,p, q) introduced in 0. 
For g = 2 it is not difficult to perform direct numerical transformation to 
construct the required operator. For g > 2, however, such an approach 
may appear to be complicated because the transforming operator is of huge 
dimension. Fortunately, it is possible to obtain simple explicit expression for 
A~^{n,p, q)Hp^q. We start with more general task of calculation of A{n,p, g) 
action on products of the creation-annihilation operators. 
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Operator A{n,p,q) was defined in |^ via its matrix representation with 
respect to specially selected basis closely related to the basis of determinant 
generators: 



(P\ (q-s) /p+\KnK'\-q~-l\ 

\q~s) \ q ) K'CN\{IUJ) V|A'nA''|/ 

(39) 

where e^j^ = | J U fsT) ( J U K| and / n J = 0, s = |/| = |J|. 

Therefore, as the first step, it is necessary to expand the restrictions of 
products of the creation-annihilation operators on g-electron sector of the 
Fock space via the basis operators e^. Using technique of manipulations 
with phase prefactors developed in 0, it is easy to derive the following 
formulas: 

KcN\{i} 



a]aj 



(9-1) 
K(lN\{i,j} 

a\a]ajai = Yl ^Ku{i,j}^ ^ + (40c) 



K^N\{i.,j} 
(9-2) 



a\a\aiai = ^ ^aS?}' ^ 7^ j ^ ^, (40(i) 

KcN\{i,j,l} 



(9-2) 

ala]aiak = ^ (-l)l^^'>"^{'>l+l^'^"^<^>lejj'^'^^'=''\ z ^ j k 1. (40e) 

KcN\{i,j,k,l} 

The second step is straightforward. One should apply the operator A{n, p, q) 
to the right-hand sides of Eqs.(40a)-(40e). Simple but somewhat tiresome 
set-theoretical and combinatorial manipulations with extensive use of two 
classic relations (see, e.g.,p[) 

V m I \k j \r — m 
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and 



lead to the following equalities: 

A{n,p,q)alaj = 



r + k — 1 
k 



P t , X 
a- a, + Oj, 



n — p n — p 



A{n,p,q)ala]aiak 



I j- P f r- j- X t 

-i-^a^aMiak I Oika-ai + ojia-ak 

(2) g V ^ 



(41a) 



(416) 



where Iq is the g-electron identity operator. Now, taking into account the 
equality A~^{n,p,q) = A{n,n — p,q) (see [|]), we can write down the exphcit 
expression for the transformed reduced Hamiltonian: 



1 (n-p) 
A-\n,p, q)Hp^q = - [Tr{h + F^)] - ^ J](z|F^|j)aJa,- 
^ U i,j 

where the Fock operator associated with some MSO index set X C iV is 
defined as 

Fx = h + J2i-Jk-Kk). (43) 

kex 

The first term on the right-hand side of Eq.(42) is the Hartre-Fock (HF) 
energy of the n-electron state |A^). It is pertinent to emphasize that the 
operator (42) should be applied to wavefunctions from the g-electron sector 
of the Fock space. Note also that with the aid of Eqs. (41a-41e) it is possible to 
expand any physically relevant operator via orthonormal basis of generators 

thus reducing the matrix element evaluation problem to the calculation 
of standard scalar products. 

The energy of the vector Z-cell in fixed MSO basis for g > 1 is obviously 
defined as 

E(J,,p,,(Z)) = min Tr [A-\n,p,q)Hj,^qdg{{i>zz'}z'cN)] • (44) 
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For g = 1 it is necessary to lift one-density operator to the level q — 2 and 
then define the energy of Z-cell Jn,p,i{Z) as 

E{Jn,p,i{Z)) = min Tr [A-\n,p,2)Hp^2d2{{iJzu{^}Z'}z'cN)] , (45) 

Wz\\=l 

where i is an arbitrary index from N\Z. 



Pcirametrizations and Induced Sheaves 

Arbitrary vector Z-cell Jn,p,q{Z) is a subspace of the vector space Sn,p,q of 
the dimension (^p"^) . Its image with respect to the isomorphism Sn^p^^ o s~^p g 
is a sum of vector Z\i^-cells 

Sn,p,.OS-l,^g{Jn,p,g{Z))^ Jn,p,^{Z\K). (46) 

K(ZZ 

(sec Proposition 8) being a subspace of the vector space Sn,p,^ of the same 
dimension (^^'^). From Proposition 9 it follows that different choice of x- 
electron functions of the type of Eq.(29) can be used to parametrize different 
subsets of the vector Z-cell Jn,p,q{Z). General scheme of such parametriza- 
tions may be described as follows: 

1. Select a certain family of x-electron functions 

(-) 

i^z\K- E 4''^l^) (47) 

SCZ\K 

with a reasonable number of free parameters; 

2. To each selected nonzero x-electron function put into correspondence 
simple {p, x)-sheaf generated by this function; 

3. Construct {p, x)-sheaf 

(5--) 

{^'z}z'cN^ Y{'4^Z\KZ'}Z'CN, (48) 
KCZ 

that, due to Proposition 9, corresponds to a certain simple (p, 5)-sheaf from 

Z-CeU Jn,p,q{Z). 
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For example, functions 

{-1)\R^^z\k\c^\Z\R\K), K c z\r 
0, K Z\R 

correspond to HF determinant \R) and single parameter Cr is used to ensure 
its proper normalization. 

It is pertinent to note that in the frameworks of this scheme families of 
x-electrons functions are used to perform local parametrization of subsets 
from Jn^p^q{Z). Of course, the cases of small x and q = n — p are of pri- 
mary interest. The equality q = n — p means that Z = N on the level q, 
^n,p,q ^ ^n,p a^d. Consequently, the above scheme is applied, in fact, for 
parametrization of p-electron states. Un—p>p then it is necessary to turn 
to the hole representation of wavefunctions and operators and work with co- 
sheaves ((p, g)-sheaves with p — n — p). Since n — p > p implies p > n — p, 
for co-sheaves the equality q = n — p can certainly be reached. 

We illustrate our approach with one relatively simple but somewhat for- 
mal example. Let 

\ y scz 

be some trial x-electron function expanded over ><-electron determinants with 
indices from the index set Z,\Z\ = p + q. Rewriting it in the form 

(<?--) (-) 

KCZ SCZ\K 

(q-xr) (j>) 

= E E (-i)'^"^^^^'4''^iAi^\i?), (50) 

KcZ RcZ\K 

where C^j^^ = (— 1)'^'~'^^\-^'c2\x\_Ri and replacing each /^-component in ex- 
pansion (50) by the corresponding simple sheaf, we arrive at the (p, >ir)-sheaf 
of the form of Eq. (48) depending on (^^*) free x-electron parameters. Sheaf 
thus constructed will be called {p, K)-sheaf induced by x-clcctron wavefunc- 
tion (49) and will be denoted as [ipz]p.!<- It is easy to see that arbitrary linear 
combination of induced sheaves is induced sheaf. 
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The next step is to derive the exphcit formulas for the matrix elements of 
the >f-electron Hamiltonian Hp g ,^{Z) corresponding to the energy expression 

^ (9--) (-) 

= -7^^ Y {ipz\KZ'\A~^{n,p,x)Hp^^\iljz\K'z')- (51) 

\ X ) K,K'CZ Z'CN 

This complicated and somewhat tiresome procedure is described in Appendix 
A. 

Since the inner product of {p, x) -sheaves (see Eq.(8)) is consistent with 
the inner product of the corresponding p-electron wavefunctons, it comes 
as no surprise that the eigenvalue problem for the x-electron Hamiltonian 
obtained from Eq.(51) should be solved in a certain p-electron metric. 

With our choice of >f:-electron functions ipzXK (see Eq.(50)) the coefficients 
in expansion of the corresponding p-electron wavefunction (30) are 

(<?--) (-) 

KCZ\B. SCZ\R 

and the normalization condition for x-electron coefficients cs takes the form 
{-) / (p) \ 

s,s'cz \rcz\{sus') J 

where the comparision |i?n A^l + |_Rn A^/j = |i?n Asas'K^o'^ 2) was taken 
into account (see Appendix A of f^). 
Let us introduce the metric matrix 

(p) 

[G,UZ)]ss' = E (-l)"^"^-«'l- (52) 

RCZ\{SUS') 

It is easy to show that the matrix (52) is nonnegative. However, it is not 
necessarily strictly positive. In Appendix B it is demonstrated how to derive 
the expression for the matrix elements of Gp^q^^{Z) as a sum of binomials. 

Matrices Gp^q^^{Z) are of purely combinatorial nature (see Appendix B) 
and to study their properties methods of modern enumerative combinatorics 
should most likely be used. The explicit description of kernals of these ma- 
trices is of primary interest. It seems to be an interesting combinatorial 
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problem with unexpectedly simple and elegant solution. Indeed, numerical 
experiments led us to the following two hypothesis concerning the properties 
of these matrices. 
Hypothesis 1. 

dimekerfG (Z)] = | ('^-") if ? - - - l(niod 2), 
mmcKer[Up,,,^[z.)i | ^ if g - x = O(mod 2). 

Hypothesis 2. If g = O(mod 2) then 

ker[G',,,,i(Z)] = cE(-l)^|{4)- 

sez 

If the Hypothesis 1 is true, then for 5 — x = O(mod 2) matrices Gp^q^^{Z) 
are strictly positive and, consequently, in this case the minimization of energy 
(51) reduces to the solution of the standard generalized eigenvalue problem 
in the xr-electron space equipped with p-electron metric: 

Hp,UZ)ipz = EzGp,,,^{Z)i;z, (53) 

Explicit expressions for x-electron Hamiltonian Hp^q^^{Z) both in the particle 
and hole representations are given in Appendix A. 

When trial function (49) is chosen to be one-electron one (just a molecular 
spin orbital) then it is necessary at first to lift its i^-components to a higher 
level K. In contrast to the situation described in Proposition 7, in the case 
under discussion the index set Z should stay unchanged. One of the ways to 
do it consists in construction of the functions 



iP) / 



(K) 
R 



\Z\K'\R), 



where \K'\ — q — x. With such choice of K'-components at the level x the 
normalization condition for one-electron coefficients takes the form 

(--1) 

E [<^P,9-2,--l(^\{^'i})]s5' = 1' 

i,jez scz\{i} 

S'CZ\{j} 

and, consequently, at the one-electron level we are faced with the generalized 
eigenvalue problem with respect to some contracted metric. We leave this 
aspect of theory to future study. 
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It is clear that for fixed Z the number of hnearly independent induced 
(p, -sheaves is equal to i^^J^ — dim,[^ker[Gp^q^^{Z)]. One can hardly hope 
that linear combinations of induced (p, >^)-sheaves may be used directly in 
molecular calculations. First, the total number of free parameters is too 
small. Second, the scheme of construction of induced sheaves does not have 
any robust physical idea behind it. The situation may be partially improved 
by turning to (p, x)-sheaf 

As„,p,^(/2) + [^z\p,>., (54) 

where 

Sn.pAR) = {(-l)"'"^^''l^'\^)}«cz'c^ (55) 
is the determinant sheaf corresponding to HF p-electron wavefunction \R), 
and A is the additional parameter to be optimized. It is easy to demonstrate 
that the problem of determination of optimal parameters A and cs,S C 
Z in expression (54) may be reduced to the eigenvalue problem involving 
augmented Hamiltonian and augmented metric: 

, "yf), ){')=eJ \ „ )('\, (56) 



where 



d 

ap,q,^{Z) = — {Sn,p,.{R) I [^z]p,.) , (57) 

d 



hp^q,^{Z) ^ ^^{sn,p,^{R)\A ^(n,p,^)\['ijjz]p,>.), (58) 

and is the HF energy of the determinant state \R). HF sheaf in expres- 
sion (54) may be replaced by any available sheaf with subsequent obvious 
modifications of Eqs.(56)-(58). 

Use of more general parametrizations in accordance to the scheme de- 
scribed at the beginning of this section is more complicated because it re- 
quires development of a certain strategy for systematic selection of reasonable 
number of free parameters c^^. Dependence of coefficients in expansion (47) 
on K results in appearance of normalization conditions of the following gen- 
eral type 



E 

S,S'CZ 



RCZ\{SUS') 



1. 
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If no additional restrictions on coefficients c^^^ are imposed tlien the corre- 
sponding metric matrix is of the form 

(p) 

R(ZZ\{SVJS') 

where C, is combinatorial zeta-function (see Appendix A). In fact, the sum 
on the right-hand side of Eq.(59) may be equal either to 0, or to ±1. 

As an example one can consider parametrization that is a direct general- 
ization of the HF case: 




K C Z\R, 
K ^ Z\R 



The number of free parameters here is (^^'^) (^) • 

Parametrizations with the total number of free parameters equal to the 
dimension of the set Ens{Sn,p,2{Z)) are of primary interest. On this route we 
may hope to approach the exact solution of p-electron problem by >tf-electron 
means that is the primary goal of the representability theory. 

It seems pertinent to emphasize once again that equations of the type of 
Eq.(53) and Eq.(56) are x-electron equations and in contrast to the standard 
configuration interaction method no expansions over p-electron determinants 
(states) appear in the frameworks of such an approach. 



Orbital Representation 



Molecular calculations are normally performed in orthonormal one-electron 
bases of molecular orbitals (MO). Turning from MSO to MO basis leads to 
somewhat cumbersome formulas but simplifies concrete calculations. Follow- 
ing Handy , we replace spin-orbital index set by a pair of orbital index 
sets (M, M) of a and f3 spins assuming that the spatial parts of molecular 
spin-orbitals with the same orbital index are identical. This corresponds to 
the so-called 'restricted' theories whereas general MSO basis embraces in ad- 
dition 'unrestricted' theories. With such an approach each {p + g)-element 
subset Z C N should be replaced by the pair (Zq,,Z^) C M x M with 
\Za\ + \Zp\ =p + q. 



23 



If the total number p of electrons is fixed and no other restrictions are 
imposed, then there is no essential difference between MSO and MO repre- 
sentations. For example, the expression for the metric matrix in terms of 
molecular orbitals may be easily derived from Eq.(52) and Eq.(C.l): 

Pc«,P/3 
(Pa+Pg=p) 

where 









, Xq + 


X/3 = X, 




1 




+ 


= X, 


- \Za 








1 -Pp^ 



\Za\ + \Zp\ =p + q. 

If, in addition, the projection Ms of the total spin is fixed, then the 
numbers Pa of a and pf^ of /3 electrons are also fixed for each split basis 
determinant \Ra,Rp)- Moreover, the q value for a given pair {Za,Zp) C 
M X M can be uniquely presented as a sum of its a and /? components: 

1 = \Za\Ra\ + \Z13\R13\ = qa + Qfi- 

Definition. Pair {Z^, Z^ is called a pair of index (g^, if l-^al = Pa+?a 
and \Zp\ =pp + qp. 

Thus, in the orbital representation for the fixed value of the total spin 
projection there are + 1 different types of pairs {ZajZp) and arbitrary 
(p,?)-sheaf {'^l^{Zc,Z0)}(^Zo,^Zp)cMxM involves, in general, g-germs with labels 
{Za, Zp) of different index. 

In MO basis the expressions for {p, g)-sheaves become rather complicated 
(see Appendix C). Here we discuss only MO representation of the metric 
matrix and the Hamiltonian involved in Eqs.(53) and (56). 

Simple combinatorial arguments lead to the conclusion that {KaiKp) 
component of arbitrary x-electron wavefunction '4)[Za,Zfj) corresponding to 
some given Mg value is expanded over determinants \So„S(}) with fixed \ Sq\ — 
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— (la ~ \Ka\ ^"^^^ I'S'/j] = = Q/B — |-f^/3|- As a result, trial x-electron 
wavefunction (see Eq.(49)) should be taken in the form 

rnin{>c,qa) 
>ta=moi;(0,>i— 

where yep = >c — k^- 

The metric matrix corresponding to such choice of trial x-electron wave- 
function is 

where Pa = i(p + 2M5) and = |(p - 2M5). 

In contrast to general metric matrix (60) its Mg components (62) turn 
out to be always degenerate. 

Ms component of x-electron Hamiltonian Hp^q^^{Z) in the orbital rep- 
resentation is easily obtained from Eq.(A.12) with the aid of the relation 
(C.l): 

I ^) [^p,q,^\{Sc,S0){S'^,S'^) - 

J2 (-l)l^^■>^^^"l-^l»^^^«l(^|/.bO[Gp„-l,,„+l,.„+l,.^,+l(^a)].„uOl.;^ 



+ 



(-1) ^ (#b)[G^P/3-i,g/3+i,-/3+i,-;,+i(^/3)]s5U{i}s;,u{i} 



^ \.Gpoi,qa,'<Oi,'<a {Zo)] 
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+- 



k,li=Za\Sa 
i,3eZa\S'a 



^_ly{i}^A{i}|+|W^A{fc}|+|{fc,OnAs„|+|{^,J■}^As^||^^^.|^■^^j 



Pc-2,qc+2,>cc+2y^+2 



{k,i}S'„u{i,j} 



jSZg\S'^ 



^[Gp„-l,q,,+l,^„+l,^'„+l{Za)]s^U{k}S'„U{i}[^P0-'i-,q0+i,>'0+i 



4 E (-1) 



{k}nAsg \+\{l}nAs^ l+IWnAg, l+IOlnA^, | 



kezg\Sg 
iezfi\s'g 

jeZa\S'a 



><[^Pa-l,?a+l,>'a + l,< + l('^a)]5„u{0S^,U{j}[^P/3-l,?/3 + l,X/3+l,>'5+l('^/3)]5^U{fe}5'U{i} 



41 E (- 



^yO}nA{ij|+|{i}nA{fcj|+|{fc,/}nAs^|+|{ij}nAs/^|^^^|_.^^ 



x[G'p5-2,g;3+2,.^^+2,..^+2(-^^)]5^U{ifc,05'U{i,i} 



The analogous expression for this Hamiltonian in the hole representation 
may easily be obtained from Eq.(A.16). 
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Numerical Examples 



To test numerous complicated formulas obtained in the preceding sections 
and to illustrate our approach, we performed several ground state calcula- 
tions for small atomic and molecular systems using Eqs.(56)-(58). All HF 
calculations to generate MO basis and list of transformed molecular integrals 



were carried out with the aid of the GAMES S program |10 



Modification of the Davidson diagonalization routine |IT[ to handle eigen- 
value problem in arbitrary non-degenerate metric (generalized eigenvalue 
problem) is not complicated and straightforward. We had certain doubts 
about possibility to use this routine for the case of degenerate metrics. To 
our surprise, after incorporation in Gram-Schmidt orthogonalization routine 
outflow for kernal vectors of the current metric we obtained stable work- 
ing diagonalization procedure. It comes as no surprise that its convergence 
properties are worse than that of the original Davidson method for the stan- 
dard Euclidean metric but it normally manifests itself only in the number of 
iterations that are necessary to reach the solution. 

In Table I the results of calculations of few simple systems in ST0-3G 
basis are presented. For B, H2O, NH2, and NH3 the particle representation 
was used {p > 2m — p) whereas for Be, and LiH transformation to the hole 
representation was required {p < 2m — p). In all calculations q value was 
taken equal to 2m — p and, consequently, only the case (Z^, Zp) = (M, M) 
was under consideration. 

In the basis of two-electron induced sheaves {k = 2) only small percent 
of correlation energy is accounted (see Table II). This percent normally rises 
when K increases, and FCI limit corresponds to >f = g in full accordance 
with general theory. It is interesing to note that for boron atom the energy 
corresponding to >f = 4 is greater than that for k = 3. It is most probably 
connected with high degeneracy of the metric matrix for k = 4 (its kernal 
dimension is equal to 35 for x = 3 and 75 for x = 4). As a result, in small 
basis used, the actual number of free parameters in energy turns out to be 
essentially less for x = 4 in compare with the case k = 3. 

We did not try to calculate more extensive systems or use better bases 
because the parametric induced sheaves even if they are centered at the HF 
origin (see Eq.(54)) will hardly become a working tool of the computational 
quantum chemistry. In spite of mathematical beauty of metric matrices 
given by Eq.(B.6) and (B.7) the induced sheaves have no reasonable physical 
interpretation and for small }< include too small number of parameters to be 
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adjusted. We used them just to illustrate the principal possibility to solve p- 
electron problems by diagonalization of x-electron Hamiltonians in specially 
selected p-electron metrics. 



Conclusion 

We have presented a new way to approach many electron problems us- 
ing >tf-electron means {k = 2,3,...). With the aid of special parametriza- 
tions it is possible to reduce p-electron optimization problem to generalized 
eigenvalue problem for Hamiltonians in spaces corresponding to numbers of 
electrons less than p. These Hamiltonians and related metric matrices are 
uniquely determined by the parametrization selected. Systematic study of 
different parametrizations and their usefuUness for concrete calculations is 
still an opened field for investigation. We have suggested two relatively sim- 
ple parametrization schemes and studied in detail one of these schemes based 
on the notion of the induced sheaves. In spite of the fact that the induced 
sheaves are introduced is a very formal way and do not have any physical idea 
behind their definition, they may be used to analyze x-electron contributions 
(x = 2, 3, . . . , g) to correlation energy. 

There exists an opinion shared by a number of notable scientists that the 
romantic period of quantum chemistry came to its end with the appearance 
of universal computer programs such as GAUSSIAN, MOLCAS, GAMESS, 
etc. As we hope, the present paper demonstrates that it is fortunately not 
yet true. 



Appendix A 

Here we describe a general scheme of evaluation of matrix elements of 
products of the creation-annihilation operators involved in the energy ex- 
pressions of the type of Eq.(51). 

Let ipz be x-electron trial wave function expanded over determinants with 
indices from MSO index set Z,\Z\ = p + q. Rewriting this function as a sum 
of its X-components (see Eq.(50)) 



(<?--) 




KCZ 
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and replacing each component in this sum by the corresponding simple {p, k)- 
sheaf {iJzXKZ'jzcN with germs 

(p) 

RcZ'n{z\K) 

where C'jf"'' = (— l)''^'~'^^'^^'cz\A'\i?, we come to the problem of evaluation of 
matrix elements ^ 

{i^z\a\^z)G = 

{q->t) (p) 
K,K'CZ Z'CN Rcz'n(z\K) 

R'CZ'n(Z\K') 

(Al) 

where a stands for some relevant product of the creation-annihilation oper- 
ators and subscript G means that the inner product is taken in the metric 
determined by the matrix Gp^g^^{Z) (see Eq.(52)). 

Changing the order of summations on the right-hand side of Eq.(A.l) and 
using technque of manipulations with phase prefactors developed in f^, we 
can recast Eq.(A.l) as 

. (-) f (p) 

{^z\a\^^z)c = E 'scs'l E 

V X / S,S'CZ RCZ\S 
\. R'CZ\S' 

Z'DRUR' J 

There are five cases to be analyzed (see Eqs.(40a)-(40e)). 
Case 1: a = a|aj. 

In this case the inner sum on the right-hand size of Eq.(A.2) is equal to 

(1 - ({i},R) E = ~ C{i},R) ( ^ _ ^ 1 > 

Z'dR \ ^ / 

where 

^ f 1 if / C i?, 
^^'^ ~ \ if / ^ i? 
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is the well-known in combinatorics zeta-function of partially ordered by in- 
clusion set of all subsets of N . By an abuse of notation, we agree to use 



Ci,_R instead of C{i},R for one-element subsets of N. In this case this function 
is identical to the characteristic function of subset R accepting the value 1 if 
i e R and ifi ^ R. 

If i e S U S' or i e N\Z then the sum over R C Z\{S U S') is just 
proportional to the matrix element [Gp^q^^]ss'- If ^ G Z\{S U S') then this 
sum is proportional to the matrix element [Gp^q-i^^{Z\{i}]ss'- The final 
expression for the desired matrix element is obtained with the aid of the 
recurrence relation (B.9): 

s,S'ez J 
Case 2: a = a|aj, i 7^ j. 

The sum over Z' on the right-hand side of Eq.(A.2) is non-zero only for 
R = Ri U {j} and R' = RiU {i}. Substitution of expansion (40b) for the 
operator aja^ in matrix elements {Z'\R\alaj\Z'\R') makes it easy to get the 
following expression for the sum under consideration: 

^_-^yi?inA{,,,;|+|{i}nA{,}|+|{j}nA{,j| j^^ - P - 1 

Since \{i} fl A{j}| + |{j} H A{,;}| = l(mod 2) for i ^ j, the required matrix 
element takes the form 

r>.\ (-) 

{^z\alaj\^z)G = —7Tr=^ J2 ^*s^s' 
^ 1 ^ s,s'ez 

(p) 

X ^ (_x)l(^iU{j})nAsl+l-RiUWnAs,|+l^?inA{,_j.}|_ 

RlU{j}CZ\S 
RlU{i}CZ\S' 
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Taking into account that Ri\J{j} C Z\S implies j ^ 5", and Ri[J{i} C Z\S' 
implies i ^ S', we can write down the final expression: 

V 1 / s,S'ez 

xCj,Z\sCi,Z\S' [G'p-l,q+l,^+l(^)]5u{j}5'U{j} ■ (^■'^) 

Case 3: a = aja^aj-aj, i ^ j. 

Substitution of expansion (40c) in matrix element {Z'\R\alajajai\Z'\R') 
immediately shows that non-zero contributions correspond to the case R — R' 
and R ^ 

\ 2 ) S,S'eZ Rcz\isus') 

Thorough analysis of mutual layout of indices i,j and relevant subsets 
of the MSO index set with subsequent use of the recurrence relations (B.9)- 
(B.IO) leads to the following somewhat cumbersome expression: 



{-/) 



{ipz\ala\ajai\ip. 



Z G 



{ipz\Gp,g,^{Z)\ipz) 



- E ^S^^' 1 i~^y^^^'^^^^^'^(hZ\iSUS')[Gp-l,q+l,^+l{Z)]su{i}S'U{i} 

s,S'ez I 

+ (-l)'^''^"^^"^''0,Z\(5u5')[Gp-l,,+l,.+l(^)]5uO}5'uO-} 
j}nAsAS' I (^ijy^z\{SUS') [Gp-2,q+2,^+2 (^) ] SU{i,j}S'U{i,j} 

Case 4: a = aial^aiai. i ^ j ^ I. 



(A5) 



Substitution of expansion (40d) in matrix element {Z'\R\a\o}jaiai\Z'\R') 
readily leads to the conclusion that this matrix element is non-zero if = 
Ri U {I}, R' — RiL) {j}, and Ri ^ i. The final expression is 

{ijjz\4a]aiai\ilJz)G = E c*sCs'Ci,z\sCj,z\s' 

\ 2 ) s,S'ez 
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X 



(A6) 

Case 5: a = ajajaiOfe, i ^ j ^ k ^ I. 

Matrix element {Z'\R\a\a\aiak\Z'\R') is nonzero only \l R — Ril} {/c, Z}, 



R' ^RiD {i,j}.We have 



{^z\al(^]aiak\^z)c 
(-) 



.l)IO}nA„|+,WnA,,j| 0) ^ 4csJ(-l 

V 2 j 5,5'GZ I 



y{fe,OnAs|+|{jj}nAs,| 



><Cfe,Z\50,Z\5Ci,Z\5'0,Z\5'[Gp-2,q+2,><+2(-Z')]5u{fc,05'U{ij}| (^-7) 

Now it is easy to derive the expression for the matrix elements of general 
one- and two-electron operators: 

n 

i) = J2mj)alaj, (A8) 

n 

0= {ij\kl)ala]aiak. {A.9) 

i,j,k,l=l 

We have 



(n-p\ S - 



(-) 



S,S'eZ jez\s 
iez\s' 



(AlO) 



32 



and 



\ 2 ) [ ij=l 



(-) 

-2 Yl "^s^s' Yl 



5,S'gZ fcGZ\s 

3£Z\S' 



n 

X 



,9+l,^+l(^)]5U{fc}S' 



U{i} 



(-) 



S,S'£Z k,iez\s 
i,jez\s' 



x{ij\kl)[Gp^2,q+2,>.+2{Z)\syj{k,i}s'yj{i,j} I (^-11) 
The Hamiltonian associated with the energy expression (51) is 



3&Z\S 

iez\s' 



y<{i\h\j)[Gp^i^g+i^^+i{Z)]su{j}S'u{i} 
+ 1 ^ (^_iy{j}nA^,y\+\{l}nAi„y\+\{k,l}nAs\+\{iJ}nAg,\ 
2 k,iez\s 

i,j€Z\S' 

■x{ii\kl)[Gp^2,q+2,^+2{Z)]svj{k,l}S>Vj{i,j} (A. 12) 

As has aheady been mentioned, for p < n — p it is reasonable to turn to 
sheaves in hole representation (co-sheaves). Reduced hole Hamiltonian (see 
Appendix B of |jl|) after its transformation with the aid of A~^{n,p, h) = 
A{n,p, x) may be written as 



A-i(n,p,,.)i/|_ = -^5^(z|%>l 
\i) i,i=i 
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+ 2M ^ {tj\kl)ala]aiak, {A.13) 

i,j,k,i=i 



where 



{i\h\j) = {-iy+Hj\h\i), (A14) 

^[j\ki) = {-iy+^+'+'{ki\tj), (A15) 

and p = n — p. Using Eqs.(A.10)-(A.ll) it is easy to derive tlie expression for 
the Hamiltonian corresponding to the energy (51) in the hole representation: 

[h;,,AZ)\ss' = l[Tr(h + Fj,)][Gp,UZ)]ss' 



jez\s 
iez\s' 



+ 1 (_l)|{i}nA{i}|+|WnA{fc}|+|{fc,0nAs|+|{M}nA5 



2 

k,iez\s 
i,jez\s' 



X (U I ^0 [Gp-2,q+2,><+2 {Z)\su{k,l}S'U{i,j} , 16) 

where \Z\ = p + q and is the Fock operator defined by Eq.(43). 



Appendix B 

With each subset Z <Z N , \Z\ — p + q let us associate the hnear mapping 

Q q 

9vAZ) : 0^n,.(^) - 0-^n,.(^), (i^.l) 

x=\ >c=\ 

defined by its matrix representation 

[G,,,,.AZ)\ss' - E (-1)1^"^-^'! {B.2) 

RCZ\{SUS') 

with respect to the determinant basis set. Here \S\ — x, and = x'. 
Diagonal (with respect to k) blocks of this matrix are identical to the p- 
electron metric matrices given by Eq.(52). We agree to omit repeating index 
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>c in writing diagonal (with respect to >c) matrix elements. Off-diagonal 
(>f 7^ blocks appear in going from molecular spin-orbitals to molecular 
orbitals. 

Let 5, and S' be k- and >f'-element subsets of the index set Z C A^, = 
p + respectively. Let us suppose that 

Sl\S' = Si < S2 < . . . Si, 

z\{s ns') = zi < Z2 < ... < Zk, 

where / = |^A^'|, and k = \Z\{S n S')\ =p + q- \S n S'\. 
The index set Z' = Z\{S H S') may be presented as 



z' = [j{s,}u{[s, + i,s,+i-i]nz'), 



(E.3) 



j=0 



where Sq = Zi — 1, = z^ + l, and the symbol [s, s'] stands for the interval 
in the set = {1, 2, . . . , n}. It is easy to show that / = h— >c'{mod 2). Cases 
of even and odd I require separate analysis, since the representation of the 
set As AS' as a union of intervals depends on the parity of 1 (see Eq.(A.5) 
from II): 



2 



U [^2i-i + 1, S2i\, if / is even, 



i=l 
[|1 



(5.4) 



U [s2i + 1, S2i+i], if / is odd. 



V> i=Q 



Using decompositions (B.3)-(B.4), it is possible to express the sum on the 
right-hand side of Eq.(B.2) as a sum of binomials. We have 



[Gp^q^:it^:^'{Z)]gg, 



E nC:), if /is even. 



i + 1 

( E ri=P) 
i = l 



j2 (-i)S'"^' n C:), ifHsodd, 

'■i.'-2 '"i+i>0 i=l 

i+1 



{B.5) 
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where 

Wi^\[si^, + i,Si-i]n{z\{sns'))\. 

i+i 

By definition, Wi > and '^Wi = \Z\{S U S')\. 

i=l 

Performing summations on tlie riglit-liand side of Eq.(A.5) first over in- 
dices Tj witli odd i and tfien over indices witli even i for the case of even I, 
and in the backward succession for odd I, we arrive to the following not very 
complicated expression 

[G,,,,^AZ)]ss' = (-If E(-1)'^(S'^'^^') (S iB.6) 

r>o \ p — r / \ r / 

In Eqs.(B.5)-(B.6) the symbol [x] stands for the function extracting integral 
part of X. 

If X = x' and S^S' then Z = and only wi = |[so + 1, Si - 1] n {Z\S)\ = 
p + q — X differs from zero which means that [Gp^q^^{Z)]gg — {^~^'^~'^) in full 
accordance with Eq.(B.2). If x = q then for S S' there is no room for 
j9-clectron subsets in the set Z\{S U S') and in this case the only non-zero 
block is the diagonal one {x = q) coinciding with the g-electron identity 
matrix. 

II Z — {1,2, . . . ,p -\- q} then the sums of interval lengths involved in the 
right-hand side of Eq.(A.6) may be presented in a more transparent form. 
Indeed, in this case Sq — 0, — p + q + 1, and for even / 

^ ' . / 

1=1 j=i 

i=0 i=l 

where Iq and li are the numbers of elements from S Ci S' belonging to 'even' 
intervals [s2i-i + 1, S2i — 1] and 'odd' intervals [s2i + 1, S2i+i — 1], respectively. 
For odd I we have 

W2i = p + g + 1 + E(-l)*Si - - lo, 

1=1 i=l 
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2 I 
i=0 1=1 

The final expression for tlie matrix elements under consideration may be 
written as 

E(-ir( ^ ), if /is even, 

E(-i)^r^'^4(-^)'^---^ if / is odd. 

L r>0 P~'^ 

(5.7) 

Directly from the definition (B.2) various recurrence relations involving 
matrix elements of metric matrices may be derived. Here we confine ourselves 
to the relations that are necessary for evaluation of relevant matrix elements. 
Let us suppose that subsets S,S' G Z are fixed and -f^d-f^l = k) is some 
subset of the set Z\{S U S'). Then 

k (fci) 

ki=OKiCK 

{B.8) 

In two most important particular cases K — {i} and K — general 
relation (B.8) reduces to 

+ (_l)IWnA,,., [Gp-i,,+i,.+i(Z)],^^^^^,^^,^ , {B.9) 

[Gp,q,:>,{Z)]gg, = [Gp,q~2,^{Z\{i, 

+(_l)l»nA,.„, [G,_,,,.^,{Z\{j})]g^^^,,^^^ 
+ (_l)lb}nA,,s, [G,_i,,,.+i(Z\{z})]^^^^.j^,^^^.^ 
+ (_l)|{..}nA,,,,, [G,_2,,+2,.+2(^)]5u{,,}5'u{,,-} ■ (^-10) 
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Appendix C 

From Eq.(6) it follows that in expansion of g-germs phase prefactors of 
the type (— are involved. In Handy orbital representation R goes to 
{Ra, Rp) and Z goes to (Zq,, Zp). Let us introduce the sets Rp = m + Rp and 
Zj3 = m + Zfs with shifted by m = \M\ elements. We have R = RaARp and 
Z = Za^Zis. Using technique of manipulation with sign prefactors developed 
in @, we obtain 

|i?nA^| = \R^nAzJ + \RanA2^\ + \RpnAzJ + \Rf3nA2^\ (mod 2). 

From Eq.(A.5) of it follows that {R/s n AzJ = and n A^^| = 
|i?Q,||Z^|(mod 2). As a result, 

\RnAz\ = \R^nAzJ + \Rf3nAz,\ + \Ra\\Zf,\{mod2). (CI) 

Certain assymetry of this relation is connected with the use of Handy split 
determinant representation with a indices always going first, 
q-electron function 

(Pc) (P,fl) 

generates simple {p, g)-sheaf with germs 

i^{Z^,Zi3)iZ^,Z'^) = 

R'^cZ'^nZc. R'^cz'^nz^ 

(as) 

From Eq.(C.2) it follows that in the orbital representation with fixed Ms 
value only g-electron functions ip^z^Zp) from the vector space Tn,q{Za^Z^ 
that are linear combinations of determinants \Sa-,Sfi) with \So\ = Qa and 
\Si3\ = {qa = \Za\ — Pa, 1/3 = — Pjs) are of actual interest. 

In the orbital representation definitions (17) and (18) should be replaced 

by 

/p(^(z.,z,)) = {{Ra, Rp) C {Za, Zp) : C^n^^n^) ^ O} ; (^4) 
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Ip{Azo.,z,)) = {Ra C : Qr^,^,) ^ for some Rf, C ; {C.5) 
I^i^z.^z,)) = [R/i C Zp : q^,,^^) ^ for some i?„ C ; (C.6) 
^.(^(z.,z,)) = U (a = «,/?). (a7) 

and each simple (p, g)-sheaf <^ '0(Za,Z3)(z' ,z>) \ is generated by any 

l " " J (Z'„,Z'g)cMxM 

of its 

/2m- - \Lp{ip^Zo.,z,))\ \ 

g-electron germs (see proposition 3). 
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TABLE I 



Results of ground state test calculations: Total energies. 









Energy, 


a.u. 






Species 


HF 


>f = 2 


>c — 3 


>f = 4 


}< — 5 


X — 6 


Be{'S) 


-14.351880 


-14.366729 


-14.373701 - 


14.403655 






-24.148989 


-24.152540 


-24.169145 - 


24.161980 


-24.189265 




LiH{^^+) 


-7.862024 


-7.864463 


-7.868083 


-7.882402 






NH2CB2) 


-54.834567 


-54.835048 


-54.839962 - 


54.850959 


-54.882707 






-55.447595 


-55.447969 


-55.449684 - 


55.457770 


-55.468587 


-55.509590 


H20CA^) 


-74.962992 


-74.965180 


-74.976319 - 


75.012500 







TABLE II 

Results of ground state test calculations: Percent of the correlation 
energy accounted. 





[Ehf 


— E^]/[Ehf — 


Epci 


X 100% 


Species 




K = 3 




K = 


5 xr = 6 


BeCS) 
B{^P) 


28.7 


42.1 


100.0 






8.8 


50.0 


32.2 


100.0 


LiH{^i:+) 


11.9 


29.7 


100.0 






NH2CB2) 


1.0 


11.2 


34.0 


100.0 




0.6 


3.4 


16.4 


33.9 100 


H2OCA,) 


4.4 


26.9 


100 
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